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Abstract. In this work we study the mathematical aspects of the develop- 
ment in the continuum thermodynamics known as the " Entropy Principle" . 
It started with the pioneering works of B.Coleman, W.Noll and I. Muller in 
60th of XX cent, and got its further development mostly in the works of G. 
Boillat, I-Shis Liu and T.Ruggeri. "Entropy Principle" combines in itself the 
structural requirement on the form of balance laws of the thermodynamical 
system (denote such system (C)) and on the entropy balance law with the con- 
vexity condition of the entropy density. First of these requirements has pure 
mathematical form defining so called "supplementary balance laws" (shortly 
SBL) associated with the original balance system. Vector space of SBL can be 
considered as a kind of natural "closure" of the original balance system. This 
space includes the original balance laws, the entropy balance, the balance laws 
corresponding to the symmetries of the balance system and some other bal- 
ance equations. We consider the case of Rational Extended Thermodynamics 
where densities, fluxes and sources of the balance equations do not depend on 
the derivatives of physical fields y 1 . We present the basic structures of RET: 
Lagrangc-Liu equations, "main fields" , and dual formulation of the balance sys- 
tem. We obtain and start studying the defining system of equations for the 
density h° of a supplementary balance law. This overdetermined linear system 
of PDE of second order determines all the densities h° and with them, due to 
the formalism of RET, the fluxes and sources of SBL. Solvability conditions of 
defining system delivers the constitutive restrictions on the balance equations 
of the original balance system. We illustrate our results by some simple exam- 
ples of balance system and by describing all the supplementary balance laws 
and the constitutive restrictions for the Cattaneo heat propagation system. 
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1. Introduction. 

With any system V of partial differential equations for the fields {y 8 (x M ), i = 
l,...,m} defined in a domain D C X of the base (typically physical or mate- 
rial space-time) manifold X there are associated several classes of objects carrying 
important information on this system and its solutions. 

First class consists of symmetry groups of different type and levels of generality 
(see [23]) and the corresponding Lie algebras of infinitesimal symmetries. Principal 
characteristic property of these symmetry groups is that the transformations from 
these groups acts on the solutions of system V transforming them into another 
solutions. 

Second class consists of the conservation or more generally, balance laws of dif- 
ferent types satisfied by all the solutions of system V . 

Relations between these two classes of objects associated with the system P are 
well know in the Variational Calculus. They are are formulated as the theorems 
of Noether type of different level of generality ([23]). The most well known for 
the Lagrangian formalism is the original Noether theorem with the complement of 
Bessel-Hagen (see discussion in 23 )), the most general result in this case is the 
result obtained by P.Olver ([23]) and, independently, by A.Vinogradov, [36]). These 
results establish the bijective correspondence between the equivalence classes of the 
conservation laws for the Euler-Lagrange system of equations and the equivalence 
classes of variational symmetries of the action functional (see [53], Thm. 5.42). 

On the other hand, in the Continuum Thermodynamics, the II law postulate 
the existence of an additional (to the basic balance system) balance law - entropy 
law = £ such that the part of the source term £ corresponding to the internal 
entropy production in the system is nonnegative. 

In the mid-60th of XX cent., it was suggested (B. Coleman and W. Noll, [5]) 
that the following statements forming the "Entropy Principle" are true for any 
thermodynamical continuum system of RET type: 

(1) Entropy balance 

(ps) it + S A xA =a (1.1) 

is satisfied by all the solutions of the basic balance system () of the ther- 
modynamical system. 

(2) Entropy density h°(x,y) is the concave function of fields y l . 
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Notice that the first statement above has sense for any system of balance equa- 
tions. Thus, one can introduce the space of the " supplementary balance laws" for 
a given balance system if, see (). 

In the terms of the balance law () the II law of thermodynamics requires (in the 
Coleman-Noll form, (see [5])) that in the absence of the outside sources of entropy, 
the Clausius-Duhem inequality 

^.+d x A{p 8 v i + ^)^0 (1.2) 

is true for all the solutions of the balance system -jr. 

Later on, I.Muller suggested ([10]) that the nonconvective entropy flux S may 
have more general form then an( j has to be defined by the requirement of com- 
patibility with the constitutive relation of the system together with the entropy 
density ft. , this suggestion can be naturally rephrased as the fulfillment of the 
requirement (1) in the list above. 

As a result of these developments, the entropy principle crystalized as the Amend- 
ment (see 22J) to the II law of thermodynamics having sense for continuum ther- 
modynamical systems of higher order (i.e. where the fluxes Fjf and sources IL 
may depend on the derivatives of fields y % up to some order). This amendment 
requires that the constitutive relations of the basic balance system -fc and 
the entropy inequality (i.e. expression of the entropy density s and the 
entropy flux S in terms of the basic fields and their derivatives) should 
be such that any solution of the balance system would at the same 
time satisfy to the entropy inequality (1.2). 

In most works it is silently assumed that the entries in the entropy balance - 
s,S,Q are the functions of x^,y l and the derivatives to the same order as the 
entries in the system -Ar. This, in particular, guarantees the removal, from the 
considerations, of trivial balance laws of the II order, see [2*5] . 

This leaves the question about the structure of the right side £ of the entropy 
balance. In continuum mechanics and in the material science the form of the 
entropy production function is of the utmost importance ([IHET]). The entropy 
"production" splits £ = S out + S" 1 into internal [real production) and external 
(outside source) terms. After this, the II law of thermodynamics crystallizes as the 
requirement that £ m ^ 0. 

This leads naturally to the problem of description of all the balance laws sat- 
isfying the first requirement above (supplementary balance laws) for a given 
balance system and modulo the trivial balance laws and extraction those of them 
that satisfy to the positivity condition. 

In Section 2 we introduce basic setting: base (space-time) space X, configura- 
tional bundle 7r : Y — > X and the balance systems -jr and collect the notations that 
are used in the paper. In section 3 we introduce supplementary balancer laws (SBL) 
for a given balance system, give examples of such SBL. In Section 4 we describe 
structures of affinc subbundlcs of 1-jet bundle J 1 (7r) — > Y corresponding to a bal- 
ance system -jc and a supplementary balance law for this system. Relation between 
these affine subbundles define the main fields X 1 and the Lagrange-Liu system of 
equations for supplementary balance laws. In section 5 we present, following the 
monograph of I.Muller and T.Ruggeri, [19] a brief review of basic notions of Ra- 
tional Extended Thermodynamics including transition to the main fields variables 
and the symmetrical hyperbolic form of the balance system. In section 6 we study 
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the LL-system for a case where main fields are functionally independent and get 
the geometrical form of the Boillat-Millcr-Ruggeri theorem on the bijection between 
the SBL with the functionally independent mean fields A 1 and the density functions 
h° with the non-degenerate Hessian satisfying to the defining linear system of the 
PDE of 2 degree. In Sections 7 we present some simple examples of solving defining 
system. In Sec. 8 we show that in the case of 2 scalar fields y 1 , y 2 in 1 + 1 space-time, 
defining system reduces to one 2 order linear PDE that has to be hyperbolic for the 
balance system to have SBL with the density h° whose Hessian is definite. In sec- 
tion 9 we study the defining system for ft, . We show that this system is generically 
subholonomic, find conditions on the flux fields F^(y) when this system is elliptic 
or holonomic. In section 12 we analyze the structure of supplementary balance 
laws for the Cattaneo heat propagation system. We describe all the supplementary 
balance laws, determine the form of internal energy for the case when there are 
SBL and its structure when the production of the SBL is positive - formulation of 
II law of thermodynamics in terms of the defining function A (i?) of the "entropy" 
SBL. 

2. Settings and Notations. 

2.1. Space-time manifold X. A state of material body will be described by 
the collection of the time-dependent fields i = 1, . . . , m} defined in a domain 
B C E n of the physical, material or ... space with the boundary dB. 

Product of closure B and the time axis T - X = T x B is the space-time manifold 
- the cylinder in the Newtonian or Lorentz space-time. We assume that the Pseudo- 
Riemannian metric G is defined in the space-time X. An example of such a metric is 
the Euclidian metric g — dt 2 + h in the Newtonian space-time ([IE]) or Minkowsky 
metric in the space-time of special relativity. We introduce (global) coordinates 
x^j/j, — 1,2,3 in B (possibly induced from the global coordinates in E 3 ) and the 
time variable t — x°. 

Denote by 77 the volume n-form rj = y/\G\dt A dx 1 A dx 2 A dx n corresponding to 
the metric G. 

2.2. State (configurational) bundle tt : Y — > X. Basic fields of a continuum 
thermodynamical theory y % take values in the space U C M m which we will call the 
basic state space of the system. 

Following the framework of a classical field theory (see [2 E] ) we organize these 
fields in the bundle 

ttu : Y -> X, X = R t x B, Y — X x U 

with the base X and the fiber U. 

Denote by s — s(x) sections of the configurational bundle - collection of values 
of fields x — > {y l (x)),i = 1, . . . , m}, x € D, where D C X is an open subset of the 
space-time manifold X . 

We will call variables y l - the vertical variables, differential of a function / G 
C co (Y) along the fibers of the configurational bundle will be called - the vertical 
differential: d v f = J^i §f d y\ 

We will also use the notation d^ for the total derivative by x^: 

, , df , , df 
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Denote by ir w : J 1 -> Y the 1-jet bundle of the bundle n : Y — > X (see [Tl [32]). 
Local chart (x ,1 ,y % ) in the domain W C Y defines the local chart (x' fi,y l , zi) 
in the preimage W 1 — ~Xio(W) C J 1 (it). Coordinates z'^ in the fibers of 1-jet 
bundle over Y are defined as the partial derivatives: zi(s(x)) = for a sections 
s(x) = {y J (x)}. 

A section s(x) : D — s> Y with the components y l (x) determine the section j 1 s(x) 
of the bundle J 1 (7r) — > X - 1-jet of the section s(x) - by the formula 

j 1 (x») = (x»,y i (x),zi(x) = ^(x). 

for a manifold M we denote by T(M) M - its tangent bundle, by T*(M) -)• M 

- its cotangent bundle. If it : Y — » X is a bundle, we denote by \^(tt) — > Y the 
vertical subbundle of the tangent bundle T(Y) formed by the vectors tangent to the 
fibers of it: V y = {£ £ T v (Y)\ir* y l; =0. In a local chart (x^, y l ) the fiber V v (tt) has 
the basis of vector fields . Denote by V* (n) — > Y - vertical cotangent bundle - 
dual bundle to the vertical tangent bundle V(ir). In a local chart (x 11 , y l ) the fiber 
V*(tt) has the basis of 1-forms dy 1 , i = 1, . . . , m. 

2.3. Balance Equations. Fields y % are to be determined as solutions of the field 
equations for the currents F?. Often (but not always) F® — y % are the densities of 
fields y l : 

F?,» = F Zt+ F tA = n - i=l,-...n. (2.1) 

Here IL(y, x) is the production+source terms of the component y % and F^(y, 

- the flow in the i-th equation. In RET theory these quantities are assumed to be 
functions of the fields y % and, possibly, of the points (t,x) — {x^} 6 X. Quite one 
restricts attention to the case where there -F/^IL do not depend explicitly on the 
space-time point x M . 



3. Entropy principle and the supplementary balance laws. 

Let 

(F?(s(x))), x » = Hi{a(x)), i = l,...,m, (*) 

be a balance system of order zero for the sections s(x) = {y l (x)} with the densities 
F®(x,y), fluxes F^(x,y),A = 1,2,3 and the source/production terms ilj(x, y). 
This system can be written as the system on the 1-jet bundle J 1 (it) using the total 
derivatives: 

d^Ff = Ui. 

Equations of system are obtained from the last equating by taking pullback by 
the 1-jet j 1 s(x) of a section s : x — > Y. 

A natural question generalizing the "entropy principle" of Continuum Thermo- 
dynamics (|19jl is - are there, except of the linear combinations of balance equations 
of the system jf, nontrivial (see below) balance laws depending on the same vari- 
ables x^jy 1 that are satisfied by all the solutions of the balance system *k. Thus, 
we define 

Definition 1. Let the -fa is a system of balance equations of the RET type. 
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(1) We call a balance law 

(K^(s(x))^ = Q(x,y(x)). (3.1) 

with the density K°(x,y) and flux K A (x,y)- the supplementary balance 
law for the system it if any solution s(x) of the balance system (it) is at 
the same time solution of the balance law (3.1). 

(2) A supplementary balance law (3.1) for a balance system it is called one 
of the entropy type if for all solutions s{x) of the balance system it the 
entropy production term Q(s(x)) is nonnegative. 

Example 1. An interesting subspace of the space SBCq of supplementary balance 
laws, including the linear combinations of the balance laws of the system (it), is 
determined by the following 

Proposition 1. Let a vertical vector field £ = ^ l d y i £ V(ir) is such that the 
condition Fjf d^ 1 = is fulfilled. Then the balance law 

a*d(t*F?^) = ?Tkr) {{CFt){s{x))), x . = (f n,)(.s(.T)) 
is the supplementary balance law for the balance system it. 

Proof. Follows from s*d(fif 7^) = d(s*fif %) = f (s)d(s*(Ff r^)+s*(if (^f)?? = 

(fn t )o S + o. □ 

Example 2. If a Lie group G is the (geometrical) symmetry group of the bal- 
ance system it , it determines the family of the balance laws corresponding to the 
elements of Lie algebra g, see f [26l [27] ) . 

Example 3. The entropy principle of Thermodynamics (see ([2H EH HI]; or 
Sec. below) requires that the entropy balance 

h% = E, (3.2) 

with the entropy density h°, entropy flux h A 7 A — 1,2,3, and entropy production 
plus source E = E s + E p of a given theory. This requirement place a serious 
restrictions on the form of constitutive relation C and leads to the construction of 
a dual system in terms of Lagrange ("main") fields A considered below (see Sec. 3 
or [19]). 

In addition to the existence of the entropy balance law, the second law of ther- 
modynamics requires the fulfillment of the condition E p ^ of positivity of 
the entropy production. Thus, the "Entropy principle" is the independent 
Amendment to the second balance law ([H]) that places severe restrictions 
to the constitutive relations, more specifically to the form of entropy density, en- 
tropy flux and entropy production. 

Example 4. System of Maxwell equations for the electrical field E and magnetic 
field H in the empty space M 3 

{ Et - VXB - ' (3.3) 
[B it + VxE = 

is the balance system (of RET type). Conditions 

V-E = 0, V-H = 0, (3.4) 

admitted for the electrical and magnetic fields in the absence of electrical charges 
are supplementary balance laws for the system (3.3). 
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4. Affine subbundles of the bundle J 1 (it) — > F,"main fields" A* and 

LL-equations. 

Consider a balance system (-fc) 

(Ffo S ) i , = n i ( S ),i = l,...,m, * 
defined on the bundle 7r : Y -3 A, and calculate explicitly the derivatives in it. We 

^ F? iXll {8) + F^(s)s^ = U i (s). (4.1) 

Here A = 1, . . . , n; /U = 0, . . . , n. 
This can be written in the form 



[E*M] = =**[ni-^]- (4-2) 



We may consider these relations as defining, for each point (x^^y 1 ) € Y through 
which a solution of the balance system if passes, the system of affine planes Ai(z), i — 
1, . . . , m in the fibers Jf,, y (n) of the bundle J 1 ^)- 

]T F? >yi (x, V)4 = ILi(x, y) - Ft iX , (x, y) (4.3) 

For each i affine planes z — » Ai(x,y) form, upon varying of the point (x,y), the 
affine subbundlc Y 3 (x,y) — > A^ C J* (77). Intersection W of these planes gener- 
ates the affine subbundlc Y 3 (x, y) — > W(z) C J* (77) (generically) of codimension 
to. 

Arguments presented above can be summarized in the following 

Proposition 2. A section s E T(D(s), Y) is the solution of the balance system if if 
and only if the image of the domain D(s) of section s under the mapping x -3- j 1 (x) 
in intersection with the fibers J^ x y \{^) of the bundle J 1 (77) -3 Y lays in the affine 
subbundle W( x _ y - ) of the fiber J^ x y ^), i.e. 

j x s{x) e W s{x) , Vx e £>(«). 

Previous considerations shows that 

Corollary 1. Let a balance system if is locally solvable at each point (x,y) G Y (see 
Appendix I), then the affine subbundles A4 and their intersection W are correctly 
defined affine subbundles of the bundle 77i : J 1 (77) — > Y of the constant rank. 



From now on we assume that the balance system if is nondegenerate 

in the sense of Appendix I. 

Let now a balance law (3.1) is defined in the same space Y as the balance laws 
of the system if. Then, it generated the affine subbundle of J 1 (77) with the fiber 
given by the affine hyperplane Ak 

j2k^4 = Q- k %»- ( 4 - 4 ) 

Since for any point z of a fiber j\ (77) over a point (x,y) £ Y there exists a 
section s : X -3 Y such that s(x) — (x, y) and its 1-jet passes through a given point 
z, the following algebraic formulation of the "entropy principle" is true. 



8 



SERGE PRESTON 



Proposition 3. Let the balance system it is nondegenerate and let 

(K»(s(x)) !fl =Q(x,y). 

is a balance law with the same domain D 1 <zY as the balance system if. 

The balance law (2.1) is the SBL for the system if if and only if for any point 
(x,y) E Y the intersection Wr XtV \ of the affine hyperplanes (4-3) is contained in the 
affine hyperplane A K ^ x ^ y) (44). 

In such a situation the following simple result of affine linear algebra is useful 

Lemma 1. Let E N be a vector space (over a field k) and let W be an affine 
subspace in E equal to the intersections of fibers of linear Junctionals fi : E — > k: 
W = rii/ i _1 (c i ), a E k. Let h : E — >• k be a linear functional on E and let c E k. 
Then the following statements are equivalent 

(1) W C h-^c), 

( 2 ) h = A*/i for some X E k. 

Applying this lemma tot the affine hyperplanes A; t and we get, the system 
of equations for a SBL (3.1) in the case of balance system of the RET type. 

Theorem 1. Let the balance system it is nondegenerate system of the RET type 
and let ok is a balance law with the same domain Y as the if. The balance law 
(3.1) is the supplementary balance law for the balance system if if and only if for 
any point z E Y there exist functions A 4 (z) E C°°(Y) (Lagrange multipliers, main 
fields) such that 

\Y.iX(*)$U-F?, x *) = Q(*)-K%»- 

For shortness we will call system of equations (4.5) - the LL-system for the 
balance system if. 

Using the differential by variables y 1 (vertical differential) we can write the sys- 
tem of equations 4.5 in the invariant form 

d v K» = X(z)d v F?. (4.6) 

First subsystem of LL-system (4.5) determines up to addition of arbitrary 
functions of . It is clear that together with all the functions of the form 
= + K^(x) solve the first subsystem (4.5) with the same functions X. 
Second subsystem defined the source term Q for the balance equation with the 
density-flux and the source term Q + for the balance law with the density- 
flux . Difference of these two balance laws will have the form d^K^ = i.e. 
is the trivial balance law of the second type independent on the fields y l , i.e. of 
order — 1. This prove the following 

Lemma 2. For a fixed set of main fields X 1 the system (4-5) determines the sup- 
plementary balance law a = K^n^ + uniquely, up to a trivial balance law of 
second type and order —I. 

Example 5. Consider a special case where all the main fields vanish: A* = 0, i = 
1, . . .,m. In this case equations (4.5) becomes statements that does not de- 
pend on K^j — 0, fj, — 0,1,2,3; j = l,...,n. In addition to this the 
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source/production term takes the form 

Q = K%». (4.7) 

Formulating now the balance law corresponding to these K 11 , Q we immediately see 
that Q — d^K^, i.e. that this balance law is the trivial one of the second kind 
of very special type ([25])- 

Example 6. For the Cattaneo heat propagation balance system (see (12.1) below) 
with dynamical variables $ -temperature and q - heat flux the LL system (4.5) has 
the form 

(X°e ~„+ TA X A q A = K° (K A =X A , 

UC+^W ; 14=^ ' A,B= 1.2,3, (4.8) 
where instead of the temperature, the variable $ = A(i?) was used. 

Remark 1. Correspondence between the supplementary balance laws (3.1) and 
the sets {X l (z)} of main fields is linear. This correspondence defined the mapping 
from the space of all supplementary balance laws of the balance system + to the 
space of sections X l d y i of the vertical tangent bundle V(tt) — >• Y. Vector fields with 
X 1 constant (^-horizontal with respect to a connection v, see |27j ) corresponds to 
the linear combinations of the balance laws of the system 

5. Formalism of Rational Extended Thermodynamics (RET). 

Here we describe, in a short form, the basic structure of the Rational Extended 
Thermodynamics developed by I.Muller and T.Ruggeri, [III [19]. For more complete 
presentation of the formalism of Rational Extended Thermodynamics we refer to 
the monograph |19j . Chapter 3. Here we introduce only necessary material in the 
form suited for our purposes. 

We will be using the notations and setting described in Sec. 2. In this section we 
take n = 4. 

Balance equations for the fields y' l (t,x) have the form 

F i!„ = F °t + F t A = n *> * = 1, • • • , m. bigstar 

Remark 2. In the rational Extended Thermodynamics one considers a case where 
balance equations are written for all the basic fields in the state space and only 
for them and where flows F- 1 and productions IT depend on the fields y l but not 
on their gradients or time derivatives. If the constitutive relations depend on the 
derivatives of fields y l , the basic state space is extended by adding these derivatives 
to the list of basic variables y l . As a result, the balance system above is the 
system of quasi-linear PDE of the first order. 

To close system of equations + for the fields y l (x) one has to to choose the 
constitutive equations of the body - to specify the densities F®, flows F A and 
production terms IT as the functions of x^,y l . 

5.1. II law of thermodynamics and the Entropy Principle. In addition to 
the choice of constitutive relations for closing the balance system, one has to satisfy 
the II law of thermodynamics. This law requires the fulfillment of one additional 
balance equation - entropy balance, with additional condition of nonnegativity of 
internal entropy production (see below). Entropy density h° is not considered as 
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an independent physical field. In homogeneous thermodynamics constitutive law 
determine entropy as the function of energy and other extensive variables (see 
[1]). In continuum thermodynamics entropy density and flux are supposed to be 
functions of basic dynamical fields y l and their derivatives. 

To facilitate the choice of constitutive relations and ensure the coherent phys- 
ical behavior of the balance system of continuum thermodynamics, the entropy 
principle was suggested. 

As we will see below, utilizing of the entropy principle facilitates the choice 
of constitutive relations for the original balance system plus the entropy principle 
and, in the case of RET allows to reduce this process to the choice of entropy flow 

3- form and to the choice of production 4-forms subject to the positivity condition. 

We start with the formulation of the entropy balance. 

Entropy density (per unit of volume) h° and the entropy flux h A ,A = 1,2,3 are 
assumed (in RET) to be functions of the basic state variables y l . These quantities 
satisfy to the balance law 

d(h»r)») = d t h° + d xA h A = E e * + S m , (5.1) 

where the right side of entropy balance is the sum of external entropy source 
Yf x and the internal entropy production E m . 

The II law of Thermodynamics is the requirement that the internal production 

4- form E m in the entropy balance is nonnegative 

E m ;> 0. (5.2) 

Entropy principle (introduced by B.Coleman, W.Noll and clarified by I.Muller, 

20 ) consists of two parts: 

(1) Constitutive relations in the balance system it and the components (h 11 , £) 
of the entropy balance should be such that any solution of the balance 
system -A- would also satisfy to the entropy balance (5.1). 

(2) (Convexity condition). Entropy density h (y) is a concave function of vari- 
ables y l : 

d 2 h° 

— — 77- — : ~ negative definite (5.3) 
oy l oy : > 

Remark 3. Conceptually, the entropy principle formulated above (see also [?]) is 
the mathematical Amendment to the physical II law of thermodynamics 

(as it was explicitly formulated and supported by W.Muschik and H.Ehrentraut in 
[22] ) . It is a very powerful tool to formulate the entropy balance for a given balance 
system and to impose limitations to the constitutive relations of the original balance 
system. 

On the other hand, as the reader will see in the next Chapter, this requirement 
is actually determine a natural extension of the original balance system to include 
all the balance equations (or, more physically reasonable, classes of equivalence of 
such balance equations modulo trivial ones) compatible with the original system. 

It follows from this that the range of validity of the first part of the entropy 
principle is much larger then just the entropy balance. 

Remark 4. Fulfilment of the convexity condition guaranteed the thermodynamical 
stability (see ([19]) or the Coleman-Noll Theorem in [35], Sec. 89) 
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Remark 5. * If the convexity requirement is satisfied, the symmetrical bilinear 
form 

d 2 h° 

9ijiv) = -Ww (5 - 4) 

can be considered as a Riemannian metric in the basic state space U. This is 
the Wundlender-Ruppeiner thermodynamical metric (|17j). It would be interesting 
to interpret the curvature of this metric in the context of RET. 

5.2. Convexity of entropy density and the main fields. Requirement that the 
entropy balance equation (5.1) is fulfilled for all solutions y = {y 1 } of balance 
system -jc leads to strong limitations on the form of constitutive equations. In 
order to apply the scheme below to more general systems we need to replace the 
fields y 1 with the new variables w l = F®, i = l,...,m. This is possible only is 
the change of variables y % — »■ w k is invertible. Thus, we introduce the class of RET 
systems satisfying to this condition 

Definition 2. Balance system of RET type (of order zero) is called regular in a 
domain U C Y if the vertical differentials d v F®(y)\yr n -\ form the coframe of the 
vertical tangent space V(tt) at each point y € Y . 

We restrict our considerations to the case of regular RET balance systems. Reg- 
ularity condition ensures that the densities w" 1 — F® of the balance laws it are 
(locally) functionally independent as functions of y 3 . As a result the functions 
w % can be taken as local coordinates in the fibers of the bundle 7r. This change of 
variables induces the change of frame, coframe and the corresponding change in the 
vertical differential d v = d y i <E) dy % — > d w j ® dw 3 . We will keep the same notation 
for vertical differential - it is clear from the context which variables are using. 

Notice that the conditions of the Entropy Principle, including the convexity 
condition of the entropy density (- positivity of the minus Hessian of the function 
h ) are invariant with respect to the change of variables y — > w. 

Remark 6. Introduced regularity condition guarantees that the balance system 
can be written in the normal form 

d t ... +d X A ... = .. 

VW \ F nJ V n « 

and, therefore, is nondegenerate as the evolutional system. If differentials dF®, i — 
1, . . . , m are functionally dependent on the fibers of the bundle tt, system describes 
evolution of some but not all fields w l and leaves some other fields or combination 
of fields to satisfy some stationary equation. 

The first condition of the Entropy Principle, formulated above, is equivalent to 
the following two statements: 

1. There exist functions X l (w) (called Lagrange multipliers or "main fields") on 
the space Y such that for all values of variables w % 

dh» ■ dF? 

— — - = X 3 ■ -—-^r dh^ — X 3 ■ dF? (5.5) 
aw 1 aw 1 J 

and 

2. E = A* II; ;> 0. (5.6) 
For the proof of these statements we refer to Section 4 or |19j . Ch.3. 



12 



SERGE PRESTON 



First of the equation (5.5) defines the Lagrange multipliers 

dh° 

dw* ^ 
Introduce the m-dim manifold A m with the coordinates X l ,i = l...,m and the 
trivial (dual configurational) bundle ir d :IxA->I. Then the correspondence [5T71 
can be considered as the gradient mapping from the fibers of the bundle tt : Y — > X 
to the fibers of the dual bundle ir d . To define this mapping in an invariant way 
identify the manifold A with the fiber(s) of the dual vertical bundle V*(ir) — > Y — > 
X . This mapping is then the section of the bundle V*(ir) — > Y defined by the 
function h° . 

Differentiating last relation by tip we get 

Hess(h°) = & ' — = — 

To qualify the invertibility properties of introduced gradient mapping y — > A we 
use the following result ([9],Ch.7, Sec. 1.1. 

Lemma 3. Let D C R N be an open domain in R N with the coordinates x l and 
f G C k (D), k ^ 2 is a smooth function in D. The gradient mapping <j> : D — » M. N : 
cj ) i (x) = §l. Then 

(1) The gradient mapping (f> is locally invertible if 

det(Hess(f)) ^ in the domain D. (5-8) 

(2) If the domain D is convex and the Hessian matrix Hess(f) = [f^'xi) * s 
positive (or negative) definite on D, then the gradient mapping (j> is C k ~ 1 - 
diffeomorphism of D onto D* = <j)(D). 



Proof. If the propertv l5.8l holds. then, due to the Inverse Function Theorem, locally, 
(f) is a diffeomorphism of the class C fc_1 . 

Let now D is convex and the Hessian of the function / is positive definite at 
all points of the domain D. We have to prove that cf> is one to one in D. Let 
<p(xi) = 4>(x2) for some x\,x% £ D and let x = X\ — X2- since D is convex, 
points x% + tx, ^ t ^ 1 belong to D. Then the correspondence t A(t) = 
Hess(f)(xi +tx) for t e [0, 1] defines a continuous matrix valued function with the 
positive definite A(t). We have, using euclidian scalar product in R , 

=< x, 4>(xi) - 4>(x 2 ) >= (x, [ ^-(i>{xi + tx)dt) = [ ,x,A(t)x>dt. (5.9) 

Jo at Jo 

Positivity of matrix A(t) for all t shows that x = 0, i.e. x\ = x%. □ 

Corollary 2. If the entropy density h° is a strongly convex (or concave) 
function of its arguments w l , then the change of variables \5. 7\ {w' } — > {A J } 
is globally invertible and defines the diffeomorphic mapping 

p:U^A (5.10) 

from the state space U onto the space A = p{U) C R m of values of variables 
A = {A*}. 
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5.3. Dual formulation. Using introduced change of variables, one may present 
all the quantities as the functions of dual (Lagrange) variables A 1 : 

Ft = ^"(Aj.fli = fii(A); h» = fc"(A). (5.11) 

Introduce the four-vector potential (3-form in 4-dim space-time) 

h = fr% = [A* • Ff - h"(X)]v„ = AiFi - h., (5.12) 

where last equality is the definition of 3-forms Fi.h. 

In terms of functions the relation (5.5) takes the simple form 

dh"(X) = F?d\\ (5.13) 

here summation by repeating indices is assumed. 
From the relation (5.12) it follows that 

Ff = |- = A* (5.14) 

As a result, 4n+4 constitutive functions F? and ft. M (A) can, in terms of A variables 
be derived from the 4 functions - coefficients of 3-form h. This representation 
is dual to the relations d„/i M = X 1 d v F] J ' (see (5.5)). 

Taking derivatives in the first equation (5.14) we get = ®^ x . ■ 

Taking in (5.12) fi = we see that the function h°(X) is the Legendre transfor- 
mation ([?]) of the function h°(y): 

f)h° 

h°(X) = ^-—w % - h°(w k ) = AW - h°(w k ), (5.15) 
aw 1 

where w % is considered in the last term as the functions of A* obtained by solving 
equations (5.7). 

Lemma 4. If the Hessian Hess(h a ) of the function h (w) is negatively definite, 
so is the Hessian Hess(h°) of the function h°(X). 

Proof. From dh° = X l dw l + y l dX l - h° wk dw k — w l dX l , we see that up — Using 
this we get 

d 2 h° _ dw l _ dX l x _ d 2 h° ! 
dX'dXi ~ dAi ~ ~ Wdw^ < ' 

□ 

After presenting currents F^ in the form (5.11) what is left of the requirements 
of entropy principle (provided the condition of convexity of h is fulfilled) is the 
residual inequality 

£(A) = A'lT ^ 0. (5.16) 

Two statements containing here determine the entropy production S in terms of 
the production 4- forms IT and require non- negativity of S. 

Reversing the arguments leading to the statements (5.14) and (5.16) one proves 
the following basic result of RET leading to the dual formulation of system of 
balance equations and the entropy principle (5.1) 

Theorem 2. |19) The following statements are equivalent under the condition of 
the convexity of entropy density h (w l ) as the function of fields w l : 
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(1) Entropy principle is fulfilled for the balance equations -fc and the entropy bal- 
ance equation (5.1) for given constitutive functions Ff'(y),Ili{y),h^(y),Ti(y). 

(2) Constitutive fields ' (y),h^{y),Y,{y) are obtained by the relations (5.7), (5. 14), 
(5.15) from the four-potential h(X) (formal 3-form) and the production 4~ 
forms IIi(A)ry for which the residual inequality 

a 1 il ^ o 

lis fulfilled. 



2. DEFINE in Appendix the symmetric hyperbolic systems. 



Combining balance equations it with the change of variables w l — > A-? and using 
the relation 5.14, we rewrite the balance system in the form 

dF^ dX j ~ d 2 h^ d\ j 

This is the system of first order PDE for the fields X i (w) (of the type first suggested 
by S.Godunov, see [Ml HU) with the symmetric matrix gfi^j- Lemma 4 proves 



d 2 h^ 

im ii.z,j i wiLii Lilt: s^iimicuiiu iiicilii^. ■ 

that the matrix awSw is negatively definite. This finishes the proof of the following 



lemma where symmetry of the matrices g^jg^j is obvious. 

Lemma 5. Let the function h is concave or convex and let w 1 —> A 3 is the change 
of variables w l to the main fields A 1 . In the variables X 1 , the RET balance system 
~fc has the form of symmetric hyperbolic system (by Friedrichs, see [8],) 

for the Lagrange fields A 1 . 

6. Main fields: Case of functionally independent A 1 . 

In the RET case the Lagrange-Liu system of equations (4.5) for the fluxes and 
the source of the SBL of the balance system t*t has the form (we are using variables 
w l considering system regular (see Definition 2)). 

j£i X i (w)d v F^ = d v K», & J2i = K %*> V/i = 1, . . .,n; j = 1, . . . ,m, 

(6.1) 



with d v being in this case the vertical differential in the bundle n: 

dvf(x,w) = f, w idw\ 

In this and the next section we omit the sign of variables in the arguments of 
functions silently assuming that this dependence is present. No derivatives by 
will be calculated in these sections, so this will not lead to any ambiguities. 

Fix an index fi in the first of these relations and notice that locally in Y the 
exactness of the form in the left side is equivalent to its vertical closeness i.e. to 
the following integrability condition 
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d v (^2 y{w)d v Ft) = d v X l {w) A d v F?(w) = 0. (6.2) 

i i 

Here we restrict our considerations to the case of regular RET balance systems 
(see Definition 2 above). Regularity condition ensures that the functions w l = 
Ff(y) can be taken as local coordinates in the fibers of the bundle tt. We will keep 
the same notation for vertical differential with respect to variables w l - it is clear 
from the context which variables we are using. 

As the next step, we apply Cartan Lemma, see Appendix 3, to the second form 
of the condition (6.2) for jj, = written in terms of the vertical variables w l : 

^dyX^w) A dw l = 0. 

i 

We get: d v \ l — J2j Aijdw J , A^ = Aji — > \ % wj = Aij = Aji = \ J wi for the functions 
Aij(w). Applying the mixed derivative test we finish the proof of the following 

Lemma 6. If the RET constitutive relation C is regular (functions F® are vertically 
functionally independent) then the condition (6.2) for fi = is fulfilled if and only 
if locally (and in a domain W C Y with star-shaped fibers over X , globally) there 
exists a function h° £ C 1 (W) such that 

A* = — 
dw i ' 

Applying this lemma to the condition (6.2) we find that Lagrange multipliers 
X t (x 7 y) have the form 

dh° 

for a smooth function h° e C 1 (H / ). Substituting these expressions into (6.1) wc 
find that one should have 



'K° =h° + K°{x), 

< d v K A (y) =J2 i ^\ wi=F o {y) d v F l A ,yA = 2,...,n, (6.3) 

These relations defines uniquely the production term Q(y) in (3.1). Flux terms 
K A , A — 2, . . . , n in the balance law (3.1) can be found up to the adding of an 
arbitrary function of x provided that for all A the 1- forms on the right are closed 
with respect to the vertical differential d v . Closeness condition(s) form the set of 
necessary and locally sufficient conditions for a function h° to generate a secondary 
balance law of the system Specifically, these solvability conditions have, in terms 
of variables w l = Ff(y) the form (here d v F° = dw J and F A (w) = F A (y(w))) 

f)h a f) 2 h° 
dvg-(w k - F»)Ad v F t A = g-g- (F%)dw> Ad v F A = 0, WL = 2, . . . ,n. (6.4) 

Here F A {w) = F t A {y k = y k {w)). 

As a result we have proved the following 



16 



SERGE PRESTON 



Theorem 3. Let (-fa) be a balance system of the RET type with the functionally 
independent by vertical variables y l density functions Ff . 

There is a bijection between the regular supplementary balance laws (3.1) of 
the balance system (if) defined in an open subset W C Y such that the multipliers 
X 1 = are functionally independent and the smooth functions (potentials) 
h°(w 1 , . . . , w m ) G C°°(W), modulo addition of an arbitrary function w*f(x),f G 
C°°{'k{W)), with the nondegenerate vertical Hessian Hess(h°) — d ®ig w j such that 

d 2 h° 

g^j(w k = F°)d v F°(w) A d v F A (w) = 0, A = 2,...,n. 
This bijection is given by the relations 

'K°(y) =h°(F°)(y), 

< d v K A (y) = £ 4 XdyFf, VA = 2,...,n, (6.5) 

Mv) =^ l My) + K%^ 

where Lagrange multipliers X 1 are defined by 

dh° 

xl (y) = Q^\^=F° k (y)- ( 6 - 6 ) 

Example 7. (1) It is clear that all the linear functions h°(w) = a l Wi + b satisfy 
to these conditions. They correspond to the linear combinations of balance 
laws of the system with constant coefficients. 
(2) Let now h°(w) = ^kijW l w : ' be a homogeneous quadratic function of its 
arguments (fcy el - const). Then, the condition of Theorem 3 takes the 
form 

kydw 1 A d v F A (w) = ^ {k ij F A l )dw j A dw l <s> k l3 F A t - k a F A = VA,j, I. 

In terms of the quantities D A = kijF A , last condition takes the form D A t = 
Dfj. This condition can be rewritten in the form of 1-form closeness 

d v (D?dw>) = 0=> locallv D A = ^-t, (6.7) 

for some functions q A (w). Returning to the functions F A we get simi- 
lar constitutive condition for existence of a SBL with the quadratic by w l 
density h (w) and the local representation for linear combinations of these 
functions (and for the functions themselves if matrix kij is nondegenerate) 
in terms of "potentials" q A : 

dvikijFfdvt) = => ioca "» k tl F A = ^Lj. (6.8) 

Next we will see that the situation existing for the quadratic density h° (w) 
extends to the general concave (and convex) by w function h°(w). 

Let now h° be a concave function, i.e. let the Hessian of h° be negative definite 
in its domain 

d 2 h° 
dw l dw : > 

Q7 

Then, the transformation w % — >• X° = ^jr\w k =F°(y) is the globally defined diffeo- 
morphism. Using X 1 instead of w % — F°(y) as the new variables along the fiber of 
the bundle ir : Y — > X we write integrability condition (6.2) in the form 



SUPPLEMENTARY BALANCE LAWS 



17 



dX l Ad v F A (X) = 0, (6.9) 

where F A (A) = F A (w(Xj). Using Lemma 2 we see that this condition is (locally ) 
equivalent to the existence of the functions (potentials) h A (X) defined uniquely, up 
to an addition of functions of x such that 

Plh A 

*S*W = (6-10) 

Definition 3. Let h (w) be a concave function of variables w J such that the con- 
dition (6.4) is satisfied. Then the n-form 

h = = h°ria + h A n A 

for which condition (6.4) is fulfilled is called a (local) n-potential of the balance 

system Balance law (3.1) corresponding to the function h° will be said to be 
generated by the n-potential h, [30J . 

In a case where h° is a nondegenerate at a point (a;, wq = F°(x, yo)) in the sense 
that its vertical Hessian 

Hess(h°)(x, wo) ~ 

is the nondegenerate matrix, the change of variables {w k } — > {A J = exists 
locally, in a neighborhood of the point (x, w ) and the conclusions done for the 
convex h° will follow locally. As a result we get the following 

Theorem 4. Let in a neighborhood of a point yo the differentials d v Ff{y) form a 
coframe in the fibers Y x , i.e. functions w k (y) = F®(x,y) are functionally indepen- 
dent. Then, 

(1) Locally, in a neighborhood of the point yo there is a bijection between the 
supplementary balance laws (3.1) for the balance system -^r such that the 
main fields X 1 = Q~ are functionally independent and the functions h°(w k ) 
of variables x,w k = F®(y) with the nondegenerate vertical Hessian at the 
point Wo — F®(yo) defined up to an addition of an arbitrary function of x 
and such that for the variables X 1 — ^r\w k =F°(y) that are defined in a 
neighborhood of the point yo the integrability conditions 



ff 2 h l> \ 

dvJduJj lx ' w=wo 



dX l Ad v F A (X) = Q^d v (F A dX l ) = 0, A = 1,2,..., 11, (6.11) 

where F A (X) = F A (x,y(w(X)j) , holds. 
(2) For a given supplementary balance law (3.1) the function h°(x,w) is defined 
by the condition 

h°(x,w) = K°(x,y(w)), where w k (y) = F°(x,y), 

conditions (6.4) are fulfilled, local vertical coordinates X 1 = §^\w k =F°(y) 
that are defined in a neighborhood of the point y satisfy to the integrability 
condition and there exists functions h A (x,X l ) such that 



F A (X) = F A (w(X)) = 
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thus defining the n-potential h(X) = h^(x, Xjrj^, corresponding to the balance 
law (3.1). 

(3) Vice versa, for a given a given function h°(w k ) of variables x,w k = F®(y) 
satisfying to the conditions above at the point u>o — F®(y ), the supplemen- 
tary balance is defined from the conditions: 

'K°(y) = h°(F°(x,y)), 

< d v K A {y) = Y.i XdvF t A , WL = 1, . . . ,n, (6.12) 

Mv) =Ei^i(y) + K%^ 

where Lagrange multipliers X 1 defined by X l (y) — $£i\w k =Fg(y) are func- 
tionally independent. 

Proof. Al statements except the last one were proved before the theorem. Last one 
follows directly from the Theorem 3. □ 

Remark 7. To compare results obtained in this section with these of the RET in 
Sec. 5 we have to replace by K M here, h 11 in Section 5 by here. 

7. Examples. 

Example 8. Consider the case of one field u(x, t) in the 1 + 1 space-time. Balance 
system takes the form 

u t + c(u), x — U(t,x,u). (7.1) 
System of LL-equations takes the form 

\ X - l = K ^ (72) 

Substitute A from the first equation into the second and get 

c'{u)K° u = K] u . (7.3) 

Integrating we get 

/U pU 
c'{u)K° u du + K 1 {x,t)=c'{u)K <} - / c"(u)K°du + R\x,t). 

Thus, we have proved the following 
Proposition 4. Let 

u t + c(u) !X =U(t,x,u) (7.4) 

be a balance equation for scalar function u(t,x) in (l+l)-dim space-time. Then an 
arbitrary supplementary balance law for this equation has the form 

d t K° + d x K x = Q, 

where 

(1) Density K°(t,x,u) is an arbitrary function of arguments t,x,u, 

(2) Main field (only one here) is X — K° , 

(3) Flux K 1 is given by expression 

K 1 =c'(u)K° - J c"(w)if rfu + if 1 (a;,t), 
where K(t, x) is arbitrary, 
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(4) Source/production term Q is 

Q = K° ■ II + K° t + K] x . (7.5) 

Example 9. As the second example we consider the Cattaneo heat propagation. 
Original basic fields are temperature y° — ■& and the heat flux y A = q A , A = 1, 2, 3. 
Balance equations 

(d t ( P e) + d X A(q A ) = o, 
\d t (r q B ) + d xA (s A m) = -q B - 

Here r = r(i?) ^ is the relaxation time and A(i?) is the function such that in 
the expression d xB (A({)) = A(i?)^^-, A(i?) - heat ... coefficient. We assume that 
p — const. 

For simplicity we consider here the case where r — const, more general case is 
considered below, in Sec. 12. If r — const, second equation can be written in the 
form 

d t (q B ) + d xA {5 B A T- 1 m) = -T- l q B . 

Introduce variables w° = e = pe, w A = q A . In these variables, the balance system 
(7.6) takes the form 



d t (e) + d xA (q A ) = 0, 

d t (q B ) + d xA (S B A(i,q c )) = -r- 1 q B 



(7.7) 



with the new function A(e, q c ) = t 1 A(?9). 
We have: F A = q A , F A = 5 A A. 

System of equations (6.4) takes in these notations, the form 

\{*l A w kd w jd<: - q A wj d w kd<; + (5%A)^ k d wJ d q c - (5%A) tWj d w kd q c]h° = 0, 

A,C = 1,2,3; k^ j;k,j = 0,1,2,3. (7.8) 

Take, at first, k = 0,j = B > and, correspondingly, w k — e,w B = q B . Then, 
previous system of equations takes the form of subsystem I. 

/ : [-S A d? + A- e d q Bd q A - A q Bd- e d q A]h° = 0, VA, B = 1,2, 3. (7.9) 

Consider separately two cases 

U = B, [-d? + A.^ qB A qB d,d qB ]h° = 0, 
[A ^ b, [AtdqBdgA ~ A q Bded q A]h° = 0. 

Take now the rest of cases (we remind that k ^ j and that the system (7.8) is 
antisymmetric by (kjj): k = D / B = j, w k = q D ,w J = q B . For this case, 
equations (6.14) takes the form of second subsystem of (7.8): 



II : [S^dgBdz - 5 A d q od~ e + S%A qD d q Bd qC - 5%A qB d q Dd q c]h = 0, 

A,B,C,D = 1,2,3; D ^ B. (7.11) 

Consider two cases: 

A^C: l5 A d qB dt~S A d q »dt}h° = 0, 

A = C: [6 A ,d qB de-6 A -d q Dde + A. qD d qB d qA - A qB d q Dd q A]h° = 0. ' 
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Let indices A = D, B,C arc all different. Then the first equation in (7.12) takes 
the form 

d qB d s h a =0. (7.13) 
For arbitrary B = 1,2,3 we can choose A,C such that the condition that indices 
A = D,B,C are all different holds. Thus, the condition (7.12) is fulfilled FOR 
ALL B. Then, the first equation in (7.12) is fulfilled due to this condition. 
At the same time, second equation in (7.12) takes the form 

[KqDdgBdqA - A q B 8 q D d q A ]h° =0, B^D. (7.14) 
Equations (7.10) now take the form 

f [-flf + A i£ -^ B ]g° = 0,B = l,2,3 ) 
[A^B, h~ e d q Bd q Aq a = 0. 

From the condition A^ ^ it follows that A ; g ^ 0. As a result, second of the 
equations (7.15) is equivalent to the condition 

d q Bd qA q° = 0, A^B. (7.16) 

In equation(7.14) indices A, B, D can not coincide. On the other hand, if all three 
of them are different, then (7.15) is fulfilled due to (7.16). This left the case A — 
B ^ D (and symmetrical case). In this case A ^ D and (7.15) reduces to 

A q Dd% B h° = 0, B^D. (7.17) 

As a result, equations (7.8) reduces to the following system of conditions: 

(A q nd 2 qB h° = 0, B^D, 
d q Bd qA q° = 0, A^B, 

d qB d~ e h° = 0, { ' ' 

U-^+A >e -^ fl ]g° = 0,B = l,2,3. 

From the second and third condition and from the invariance of h° under the spacial 
rotations it follows that 

h° = 1 (~e) + a\\q\\ 2 +(3 A q A , (7.19) 

where a, (3a do not depend on e and 0a is the spacial covector. 

In the first condition requirement A q o = for all D is equivalent to the require- 
ment that e q c = for all C. Thus, if e does not depend on the heat flux, then 
A q D = and the first condition in (7.18) is fulfilled. Then the last equation in 
(7.18) gives, for h° in (7.19) 

7"(e) = 2aA i? o 7(e) = 2a J Ade + ce + d. 
Thus, in the case where e „c =0, 



g 

A 



I 



2a I Ade + a\\q\\ 



+ cl + p A q A + d. (7.20) 



Now consider the case where A q D ^ 0, then d 2 B h° = for all B. Using this in 
(7.20) we see that in this case a = 0. Thus, we can conclude that 

Proposition 5. For the Cattaneo heat propagation model with t$ = 0, 
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(1) // l A c = 0, VC, then 



h° = 

(2) If to ^ 0, then 



I 



2a / Ade + a\\q\ 



+ cl + /3 A q A + d, 



h° = ci + (3 A q A + d. 



In the expression (7.20), last term corresponds to the secondary balance laws 
that do not depend on the fields (trivial), previous two terms correspond to the 
initial balance system (7.6) but expression in brackets gives the density of new 
balance law. See below, Sec. 12 for more details. 

Remark 8. Form of h° obtained here does not allow to determine the fluxes h° - 
for this one has to solve system of equations for h A . In Sec. 12 the whole LL-system 
of equations will be solved and, together with the description of all supplementary 
balance laws, corresponding constitutive limitations on the form of internal energy 
e = e(6, q) will be obtained. 

8. Example: (2+2)-RET balance systems. 

Example 10. Consider a model example of a balance system with one space vari- 
able x and two fields y 1 ,y 2 . The balance system has the form 

d t y i + d x F l 1 = U h i= 1,2. (8.1) 
In such a case the system of equations for h° reduces to one equation with 
k = l,j = 2: 

{Fl yl d yi d y2 -Fl y2 d yi d y i)h° = 0, 

or 

[-Fl y2 d 2 ylyl + (Fl yl - Fl y2 )d yl d y2 + Fl yl d 2 y2y2 ] h° = 0. (8.2) 

This is, for each point (t, x) the homogeneous linear equation of second order 
with two independent variables. Space of solutions of this equations is infinite- 
dimensional. 



Next we study the question - which of equations 8.2 have solutions h° with 
definite (positive or negative) Hessian Hess(h°) = (^ ^Q yJ ) ■ 

Near a point where type of this equation is locally constant, this equation can 
be transformed, by a change of variables, to one of three canonical types - Laplace 
equation (elliptic), degenerate (parabolic) and wave (hyperbolic). The type at a 
point yo = (x ,y ) is determined by the determinant of the matrix 

j{x,u) = ( 1 -^f, ^iV-^va. (8 . 3) 

(1) Elliptic case: Det(J) > 0. In such a case, equation (8.2) is locally iso- 
morphic to the standard Laplace equation on the plane. Its solutions are 
harmonic functions of y l . Hessian matrix of a harmonic function has zero 
trace - by the Laplace equation itself. As a result, sum of its eigenvalues 
- trace of the matrix is zero and a function h° constructed in such a way 
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can not be positive or negative definite in canonical coordinates. Yet, the 
change of variables that transformed equation (8.2) to the Laplace equa- 
tion is not tensorial and additional terms may change the sign of trace. 
Yet, let us fix a point y and make a linear change of variables that trans- 
form the quadratic form of equation (8.2) at the point y to the canonical 
(Laplace) form. In new coordinates (w 1 ,^ 2 ) at the point y we will have 
h° wlwl + h° wlwl = Tr{Hess(h a )(w i ) = 0. Now we use the fact that under 
a linear transformations of variables the Hessian of a function transforms 
tensorially and, as a result, its trace is invariant under the linear change 
of variables. As a result, in original variables Tr(Hess(h°))(y' t ) — 0. Thus, 
eigenvalues of this matrix have (nonzero) opposite signs and the Hessian 
can not be positive or negative definite. 

(2) Hyperbolic case: Det( J) < 0. Equation (8.2) by a local change of variables 
can be reduced to the standard wave equation (<9 2 2 — c 2 9 2 2 )/ = 0. General 
solution of this equation 

/ = /i(tu - cv) + f 2 (w + cv) 

has the Hessian 

tie 8 8U)-y c{ _ f „ + f „ ) f „ + f „ y 

Determinant of Hessian is equal to Det(Hess(h°) = Ac 2 f['f 2 . Thus, if 
functions j\ , f 2 have the same convexity type, both eigenvalues of Hessian 
(if they are nonzero) have the same sign and the Hessian form is definite, if 
functions /i , f 2 have the opposite convexity type , eigenvalues of Hessian (if 
they are nonzero) have the opposite sign and the Hessian form is indefinite. 

(3) Degenerate case Det( J) = 0. In such a case the equation can be reduced 
to a parabolic equation (d w — <9 2 2 )/i° = 0, or degenerate d 2 2 h° = 0. For 
a degenerate case any solution has the form h = hi(w) + h 2 (w)v and its 
Hessian 

-»•<*>= e ; r ? i). 

Thus, characteristic equation has the form 

A 2 - (h'{(w) + vh' 2 '{w))\ - h' 2 {wf = 0. 

Its roots 

A± = -\{h'{{w) + vh' 2 \w)) ± \J\{h'{{w)+vh>i{wW + h> 2 { W y 

are also of the opposite sign and, therefore, Hessian matrix of h° can not 
be positive or negative definite. Using the linear change of variables as 
for elliptic case and the fact that under such a change of variables the 
homogeneous second order equation (8.2) transforms into the homogeneous 
second order equation we get, in new coordinates w l equation of the same 
type as (8.2) such that at the chosen point y it has the form h^ 2w 2 = 0. It 
follows from this that at this point Det(Hess(h°))(w) ^ 0. By the same 
tensorial arguments as for elliptic case this will be true at the point y in y % - 
variables as well (Det(Hess(h )) is multiplied by the square of the Jacobian 
of linear transformation and its sign is not changing) and the Hessian form 
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Hess(h°)(y) can not be positive or negative definite (but can be positive 
or negative semi-definite !). 
These arguments prove the following 

Proposition 6. Let a balance system with one space variable x and two fields y , y 
has the form 

d t y l + d x Fl = Ui, i = 1,2. 
If there is a (secondary) balance equation (2.1) with positive or negative definite 
Hess(K°(y)), then the equation (8.2) is hyperbolic: 

Det(J(uj) = ~Fl y2 Fl yl - ±(Fl yl - F^) 2 < 0. 

Let now the balance system (17.1) is such that the defining equation (8.2) is 
hyperbolic, i.e Det(J(y)) < 0. Introduce notations a — —F\ y2 , 2b — (F^ yl — 
F% 2 ), c = F^ yl . Then the equation (17.2) for h° takes the form 

(ad 2 y2 2 + 2bd 2 yly2 + cd 2 yl 2 )h° = 0. 

Introduce the characteristic functions Ai 2 of this equation as solutions of the qua- 
dratic equation 

a\ 2 + 2bX + c = 0. 

Let functions ^(y 1 ,?/ 2 ) and tpiu 1 ^ 2 ) are Riemann invariants (see [5]), i,e, solutions 
of the equations 

<j) <y i - \i(f> tV 2 = 0, 
tp^ 1 — X 2 ip iV 2 = 0. 
Change of variables 

i = <t>(y\y 2 ), 
v = i>(y 1 ,y 2 ) 

reduces the equation for h to the canonical form 

h% v = 0, 

having general solution of the form 

h° = h 1 (0 + h 2 (v)- 

In variables y % (locally) general solution of last equation has the form 

h°(y i ) = h 1 (^y\y 2 )) + h 2 ^(y\y 2 )). 
Calculating Hessian of this function we get 

%?)■ <"» 

In this expression third and forth matrices are positive semi-definite -one eigenvector 
of each of them is zero, the other (A 2 + 1) is positive. If the graphs of Riemann 
invariants <ft, ip are convex, then one can construct examples of h° with positive or 
negative definite Hessian by the choice of h\,h,2 having first and second derivatives 
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of the same sign. For instance, if the Hessian of <f> is positive definite, one can take 
/12 = and h\ to be monotonically decreasing and convex function. 

Algorithmically this approach allows to choose functions hi for which the solution 
h° has the negative definite Hessian. 

Example 11. Let <j> be any harmonic function of its variables. Let x be the 
harmonically conjugate to (i.e. function (tp + ix)(y 1 +iy 2 ) is analytic. Then, due 
to the Cauchy-Riemann equations 



Hess(4>) 



_ Xy 



y ^y y 



'9T Xy 1 y 



and Det(Hess(<j>)) = —x y i y 2 — 4^i y 2 = 0. Choosing hi appropriately, one can 
guarantee the existence of a function h with the nonpositive Hessian. 

9. Integrability condition and the D-module A4c- 

Intcgrability condition (6.4) represents an overdetermined system of linear 
equations of second order for the function h°. Therefore, except of the obvious 
linear by w l solutions might exist only under the certain conditions on the current 
part of the constitutive relations. 

Consider a special case where densities F^(x, y) coincide with the basic fields y l . 
In this case w l = y % and condition (6.4) 

d 2 h° 

T dy> Ad v F l A {y) = 0, A=l,...,n 



takes the form 



dy l dyi 



F£ yk h° yiyj dy j Ady k = 0, A = 1,2, . . . ,n, (9.1) 



F -y lk pp ft ° = °' for a11 k^j,A=l,2,...,n. (9.2) 

This system has all the functions h° — J^k °k{x)y k linear by y % as trivial solu- 
tions, corresponding to the linear combination of balance equations of the system 
(~k). Generically there are no other solutions. This has its reflection in the known 
fact that the entropy principle - existence of a nontrivial solution of this system with 
the negative definite Hessian place strong restrictions to the form of the constitutive 
relations C. 

Introduce the vertical vector fields r\ A = F^ yk (x,y)d y i, then the system (9.2) 
takes the form 



(Vkdyi -r,fdy«)h° = 0, k ± j, A = 1, 2, . . . , n. (9.3) 

This system determines the cyclic "Dg-module Mc over the ring £ of smooth 
functions of y % containing the coefficients F A ([2]). Let now the components of 
the current are analytic functions. Consider their extension to the complcxification 
Y c of the manifold Y (see [5J [TS]). This allows us to extend Vg to the 2?o-module 
where 0{Y C ) is the sheaf of holomorphic functions on Y c . Choosing the admissible 
filtration of the algebra T>£ we get the graded module gr(JAc) = £[C\]/I where the 
ideal / is generated by the quadratic polynomials 

P$ (0 = Vf (00 " 4 (OCi = F k A r (x, y)CkQ ~ F* yj (x, y){ k Q . (9.4) 
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The characteristic variety Char(A4c) of £>-module M.c (more exactly, the relative 
characteristic variety with respect to the projection Y — > X (see [2], Sec. 1.6) is 
defined as the support of the £ [(^-module gr(M.c) - i.e. as the coisotropic subset 
of the (vertical) cotangent bundle V*(Y) — > Y where 

P£(£) =0,for all i,j, A. 

The set of expressions in the right side of (9.4) can be considered as the collection 
of coefficients of the 2-form 

with the coefficients being functions in Y. This wedge product is zero if and only 
if the 1-forms are proportional, i.e. if 

^(C) = ^ A (x,y)Q,A = 1, . . . ,n;i = 1, . . . ,m 

with some factors /i A (x,y). Left side of this equation can be considered as the 
result of application of linear operator L A with the matrix L A 1 — F A yi (x, y) to the 
(co) vector with components Q. Then the last equality takes the form 

Thus, a (vertical) covector (y, £) belongs to the support of the module gr(Aic) if and 
only if £ is the common (real) eigenvector of operators L A (y) for all A = 1, . . . , n. 
Generically there are no nonzero vector with this property and Char(A4c)( x ) — 
U x , x G X coincide with the zero section of the (vertical) cotangent bundle. 
For each A = 1 , . . . , n let 

P A (x,y)( V )= 1] (V-H j (x,y)) 

to be the reduced characteristic polynomial of endomorphism L A . Product here 
is taken over all different eigenvalues fj? of the endomorphism L A (x,y). It is easy 
to see that the operator P A (x, y)(L A ) annulate all the eigenvectors of L A (x,y) 
and no other vectors. To see it let E(/j,) be a root subspace corresponding to the 
eigenvalue [i E a(L A ). Then any operator L A — f^(x,y)I in the product above 
with fjp £ a(L A ),^ ^ [i leaves the subspace E([i) invariant and is invertible on 
this subspace. Since factors in this product commute between themselves the kernel 
of P A (x, y){L A ) restricted on E(/i) is the subspace of eigenvectors with eigenvalue 
/i. Therefore, condition for a vector ( e T* y (U^ y ) to belong to the characteristic 
variety of the 2?o-module Aic is equivalent to the fulfillment of the system of linear 
relations 

P A (x 1 y)(L A )( = 1 A=l,...,n. (9.5) 
As a result we have proved the following 

Theorem 5. (1) Let C be a a constitutive relation of the RET type with F® = 
y l . Then a function hP^y 1 ) generates a supplementary balance law with the 
functionally independent Lagrange multipliers (main fields) X 1 in a neigh- 
borhood of a point (xo,yo) if and only if it is a solution of the cyclic D- 
module M.q 

{rtfdy, - r,fd yk )h° = 0, k ? j, A = 1, 2, 3, 
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where r\ A = F A yk (x,y)d y i and such that the (vertical) Hessian 

d 2 h° 
dy l dyi 

is nondegenerate in a neighborhood of a point (xq, yo) 



Hess(h°)(x ,y a ) 



\x=x Q ,y=ya 



(2) Characteristic variety of the cyclic -module M.q is the union of the zero 
section of the vertical cotangent bundle V*(Y) — > Y and the set of lines gen- 
erated by the common eigenvectors of the linear operators L (y) : V* (Y) — > 
V*(Y), A = 1,2,3 defined by 

(L A (yKh = F£ y] ( k . 

Deeper study of the D-module A4c and the corresponding secondary balance 
laws for a RET system with a constitutive relation C will be presented elsewhere. 

10. Defining system for h°. 

In applications it rarely happens that F® — y l . In such a case one has to use 
solvability conditions (6.4) written in terms of variables w l : 

d 2 h° 

^-^-dwi AF* wk dw k = & {F£ wj d w ud w *-F£ wk d wi d w i)h° = 0,j^k;A= 1,2,3. 

This requires the reversal of relations w l = Ff(yi) to get j/ J = G : >(w k ) and use 
them to determine the flux components F^ (w) as functions of w l . 

It is easy to see that the statements of Theorem 5 are valid for this case. In 
particular, characteristic variety of the overdetermined system 10.1 is formed by 
the union of zero section of the dual vertical bundle V*(tt) — > Y and the common 
eigenvectors of three linear operators L A (w) : V*(Y) — >■ V*(Y), A — 1,2,3 defined 

by 

(L A (w)C)j = F£ (10.2) 
If three matrices L A - — F A wj are in generic position, they do not have common 
eigenvectors and characteristic variety of system (10.1) coincide with the zero sec- 
tion of the vertical cotangent bundle. This proves the first and, therefore, second 
statements of the following 

Proposition 7. (1) Generically, characteristic variety of the defining system 
(10.1) coincide with the zero section of the vertical cotangent bundle V* — > 
Y. 

(2) Generically, the overdetermined system (10.1) is holonomic. 

(3) Generically, in the case where constitutive relations C are analytical, the 
space of solutions of (10.1) is finite-dimensional. 

Third statement follows from the second one and the finiteness of the space of 
solutions of a holonomic system (see [TS] or [5]). 
Rewrite condition (6.2) in the form 



or, in the form 



(dh°\ 
\dw* ) 



d v {F( 



<^> d v 



Ad v (F*) = 0, A = 1,2,3, 



0, A = 1,2,3. (10.3) 
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Last equation can be written in the form 

d v [{nf ■ h°)d w 3] = o, 

using the vector fields n A = F A -di and the corresponded "twisted differential" d FA 
of the functions 

d FA ■ f = (vf ■ f)dw j . (10.4) 
It can also be presented in the form 

d v {J A d v h°) = 0, (10.5) 

where J A : V*(n) — > V*(n) is the endomorphism of the vertical cotangent bundle 
defined by F A -: 

J A : /J, k dw k -)■ / //,'///•'. 

Remark 9. As it follows from the first form (10.3) of defining system, functions 
h° = aiW % + a,a, h a £ R are solutions of defining system for all balance systems. 

Lemma 7. Twisted differential d FA : C°°(U X ) Q}{U X ) is elliptic iff Ker(F A ) = 
0. 

Proof. It is easy to see that the symbol a(d A F) : U®T*(U) -> ^}(U) has the form: 

A = {A 4 } -> F A X l dw j . 

Thus, a (constant) element A = {A*} is characteristic for d FA , at a point w iff 
i^-A* = for all j. □ 

Corollary 3. System (10.1) is elliptic if and only if intersection 

3 

p| Ker(F A ) = 0. 

A=l 

System (10.1) has the m + 1-dim space of "trivial" solutions - h° — aiw 1 + a 
where Oj, a are functions of x^ 1 only. 

Example 12. Systems of the form (10.1) may put a very diverse restrictions to 
the function ft, and to the constitutive relations in terms of Ff 1 . To illustrate this 
diversity we consider the simplest case when m = 2 so that there are only two 
variables w 1 ,w 2 and for all A, matrices F A - are constant. 

More then this, consider just one equation of the form (10.1) with the constant 
matrix Fij = Fij Thus, the equation for h '(w 1 , w 2 ) has the form 

d v {F?h° wk dw k )=0. (10.6) 

In the next Table we collect information about the from and type of solutions of 
equation (10.6) for different constant matrices F^ : 

From this table we see obvious relation between the type of equation and the 
action of 1-parameter subgroup of SL(2, R) corresponding to matrices M. 

Example 13. Let m be arbitrary and matrix M = F A j - constant and diagonal: 
M = diag(cii). Arrange diagonal terms so that M would take block-diagonal form 
such that each block Mi of size to/ x to/ corresponds to the subset of variables 
y J , j £ //, has the form a// TO , and a/ ^ a Sl s nel. Here [1, to] = l\ U . . . Ik is the 
splitting of the set of indices [1, to] into non- intersecting subsets. Then it is easy to 
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Matrix F i} - 


Equation 


Solutions , 


Fij = 


ill) 


,a^b 


Kw = 


/>° = iV) + <70/ 2 ), 


F — 


(0 a\ 
{0 Oj 




h lw = 


h° = /(»V + g{y% 


TP. . — 
r V — 


{0 1 

l-i oj 


^ v + = o 


Harmonic functions on the plane, 


F — 


(0 1 
ll o) 




/l ! 1 - /A 2 = 

^2 


h° = f{y l +y 2 ) + g(y l -y 2 ). 


Fij = 


( a 

\-b - 


-«) 


/l ll + 2^? 2 + ^ 2 =0 


Elliptic if ^ < 0, hyperbolic if £ > 



Table 1. Type of solutions of equation 10.6. 



see that each subsystem of the system (10.1), corresponding to the fixed value of 
A, splits into the equations 

"■i^Av* ~ "''P'.'// y = ( a 3 - a k)h a jk = 0. 

General solution of this system has the form 

h° = ^2h l (y1,jeli) (10.7) 
i 

where all the functions h l {y^,j 6 /;) are arbitrary functions of their variables. 

Example 14. Let M = diag(fi(y 1 )) - diagonal matrix with coefficients depending 
only on y 1 . It is easy to see that in the case where diagonal elements /, of M may 
depend on all the vertical variables, equations 10.1 have the form 

difih^dx*) = & (f t , 3 h°, + Uh%)da? A dx l = 0. 

This can be rewritten in the form 

/,,/<:: - hihPj = (fj - h)h% , i + j. (io.8) 

In the case, where coefficients /, depend on y 1 only, this system splits: 

f/i,ifc?i = (/i-/i)ft?ii, i>l,J = l, (1Q 9) 

\(fj - fi)h,ij =0, i,j > 1. 

First equation can be written in the form 

f h° 



h - fi ^ 

Integrate by y 1 and denote a antiderivative of j^zj by F l (y 1 ). We get 

ln(\h° i \) = F i {y 1 )+G i {y i ). (10.10) 

Function G l may depend on all the variables y % except y 1 . 
Last relation can be rewritten in the form 

h% = ±e Fi ^e Gi W\ (10.11) 

notice that h° — is also solution of the first equation above. As a result, we get 

fc? B «=0V)e Gi(fli) , (10.12) 
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where (^(y 1 ) are arbitrary functions of one variable y . 

Calculating now derivative by y 3 , j ^ i, 1 and substituting obtained relation into 
the second equation in 10.9 we see that 

fi U => g% = o. 

Restrict now to the case where fi ^ fj. Then for all i,j > 1 such that i ^ j, 
G 1 ^ — 0. This, finally gives us 

h° = Y J <t>\v 1 )e Gi{yi \ (10.13) 
»>i 

for arbitrary functions <jf (y 1 ) , G % (y l ) . 

This Example illustrates the fact of the following pattern; if entries of constitu- 
tive relations depend on some variable in non-symmetrical way (temperature and 
density are obvious examples), entries of supplementary balance laws will depend 
on the same variable in non-symmetrical way. In the case of Cattaneo heat propaga- 
tion (see below, Sec. 12) such variable is temperature. As a result, in the expression 
of the internal energy and in the entropy balance there appears the new constitutive 
function A°(tf). 



In practical situations it is more convenient to use variables y % because using w % 
one has to express as functions of w l and that leads to the complex expressions. 
That is why we will rewrite the compatibility equation - defining system (61.3) in 
terms of variables y l . We will use the relations 

n , , dh° dh° dy k 

l,y3 ' dw l dy k dw l 



Introduce the matrix W\ 



and assume that this matrix is non- degenerate 



(regular case). Let W 1 be the inverse matrix. Then we have in (10.1): 



^4 dh° , , ^ A dh° dy k ^ n , ; r-,4 dh° TTr , , , , 



Exterior differential of a function or forms is invariant under the change of variables. 
Thus, condition (lo.l) can be written in the form 



d v 



(V 1 k F*W?) °£ d j 



d v 



A k 



—dy l 
dy k 



d„(tf£ -d v h°) = 0. (10.14) 



Here is the vector of (l,l)-tensor fields in the vertical space, for each A = 1, 2, 3 
defines the endomorphism of the vertical cotangent bundle V* — > Y. is the 
" constitutive" quantity generalizing tensor of characteristic velocities in the theory 
of hyperbolic waves. It is defined by the constitutive relation C. 



Remark 10. In these variables linear solutions have the form h° = aiF®(y) 
a, at, a 6 R. 
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Example 15. Consider the system (10.1) for Cattaneo heat propagation balance 



system (7.6). We have here y° 



We have Det{Ff j 



w = f; 



'■(re. 



■o 



( M 

TM A,i ,q A ) 



q A ,A 



-,Q L 
T 







1,2,3. Therefore, 

,q 2 e .q 3 \ 



T 



r 



(10.15) 



Thus, the regularity condition for 



constitutive relation of Cattaneo balance system is 

T 2 (re^ - T^q A e. qA ) ^ 0. 



(10.16) 



Calculating the matrices F A j we get 



F 1 ■ — 

± 1,3 



I o 

A(i?),tf 













F 



V 



o 







F? 







\A(0),* 






A(0),* 


(10.17) 

It is easy to see that rank of each of these matrices is two, their kernels do not 
intersect: f] A Ker(F A j) — {0}. As a result, defining system for Cattaneo 
balance system is elliptic. 

More then this, square {F A j) 2 of each of this matrix is equal to times the 
projector to the two-dim plane < q A > in the vector space U =< fl,q A ,A = 
1, 2, 3 >. Derivative A^ is the heat conductivity coefficient that is nonnegative due 
to the residual entropy inequality (Muller, TD, Sec. 1.3. 2). 

In the basis eo = &a = q A eigenvalues of ■ are eo ± -i/A^e^, e^+i, e^+2 



(cyclically counted) and as a result, these matrices have no common eigenvectors. 
As a result, defining system (10.1) for Cattaneo system is holonomic. 



1\ 



oy 



11. Defining system of the 1st order. 

In this section we rewrite defining system (10.1) as the system of equations of the 
first order for the generating function h° and 3 additional functions h A , A = 1,2,3. 
Defining system (10.1) 



F A h° = F A h° 

r i,wi"',w i w k r i,w k ,w % w3 



can be written in the form 



or 



d w k (F A wJ h° i ) — d w , {F A wk h° 



This condition is locally equivalent to the original LL-system (6.1) - condition of 
existing of the potentials h A such that 

F l * wi h% i =d wj h A , A=l,2,3;j = l,...,m. (11.1) 

Last condition can be understood as the relation between the vertical differential of 
the function h° transformed by the endomorphism of the vertical cotangent bundle 
7ry : V(ir)* — > Y and the vertical differential of the potential h A : 

F A d v h a = d v h A , (11.2) 
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where T A is the endomorphism of the vertical cotangent bundle defined by 

T A {dw') = F A wk dw k = d v F A . 

Multiply the yl-th equality by the real parameter ca and add them. Then, intro- 
ducing the endomorphism F{e) = caF A of the bundle V(n)* — > Y and the function 
h(x,w 7 e) = €Ah A , we can write the last equality in the form 

T(e)d v h° = d v h(e,x,w) F{e) t . wl h° wl = h(e) tW j. (11.3) 

Here T{e)i — caF a . Vector e £ R 3 is the parameter. Last relation should be 
fulfilled for all values of e. 

Definition 4. Balance system (+) will be called C-regular if for some vector e* £ 
R 3 , the matrix F{e)i tW 3 — tAF A wj is nondegenerate: det(J 7 (e*) i . w j) ^ (and the 
corresponding endomorphism is reversible). 

If the condition of this definition is fulfilled for some e* , denote by J r ~ 1 (e*) - the 
inverse endomorphism of vertical cotangent bundle. Applying this endomorphism 
to (11.3) we get the formula for the differential of the density h°(x,w). 

d v h° =J r - 1 (e*)d v h(e*). (11.4) 

Let now we assume that the balance system is C-regular and try to reverse the 
arguments leading to the formula (11.4). Namely, let h A (w) be some functions and 
let the m x m-matrix T{e) be invertiblc for I in an open set £ C R 3 . 

Then, there are two questions: 

(1) Is the expression J 7 ~ 1 (e*)d v h(e*) the vertical differential of a function h°l 

(2) Under which condition the function h° does not depend on the parameter 
el 

The answer to the first question obviously is: there exists a function ft, such that 
for e* equality (11.4) is fulfilled iff 

d v (F-^e^dvhie*)) = 0. (11.5) 

In this equality the matrix of 1-forms acts on the 1-form. 

Assume now that (11.4) is fulfilled for some e* £ E. For correctness, the last 
result should not depend on the choice of e* and, therefore, equality (11.4) has to 
be fulfilled in the whole connected component of the set £, containing e*. 

This condition is fulfilled at least at the named connected component iff the 
derivatives by €a of the right side in (11.4) are identically zero in the whole con- 
nected component. This leads to the system of equations 

± A (F-\z ) d v h{e))=0, (11.6) 

for h A whose solvability conditions delivers the constitutive restrictions on the 
balance system Calculating derivatives here and using the fact that for a 

matrix function F(e), d^F^ 1 = —F~ 1 d f FF~ 1 , we see that the condition 

- T-H^QpiH^H^dyhie) + J- 1 (e)^^ft(e) = 0, e £ E (11.7) 
or, applying T(e) to both parts and calculating derivative in the last term, 

fA^iddvhie) = d v V A , A = 1,2,3. (11.8) 



32 



SERGE PRESTON 



has to be fulfilled. More specifically, last condition has the form 

^F^i^eBdyh 3 -d v h A = Q,^ (J" A J" _1 (e) e B - #b ® I m )d v h B = 04 = 1,2,3, 

(11.9) 

or, the form 

{J= A e B -5 A (8) FiefiF^ffidvh 3 = 0, A = 1,2,3. 
These considerations can be resumed in the form of 

Theorem 6. Let it be a RET balance system of the RET type and let (3.1) be a 
supplementary balance law such that corresponding main fields X 1 are functionally 
independent. 

(1) Then the density hP of this SBL satisfies to the system 

F* w ih%i = d wj h A , A=l,2,3;j = l,...,m (11.10) 

for some functions h A (w % ). 

(2) If for some real parameters e = {ca^A = 1,2,3} the matrix J-(e)i, w j = 
tAF A w j is non- degenerate, then 

d v h° = ^(jtfahfjt), (H-H) 

where h(e*) = tAh A - 

(3) Functions h A {w l ) generate the SBL for the balance system (if) if and only 
if they satisfy to the compatibility system 

(d v (F- 1 (e*)d v h(t*))=0, (uu) 
\(F A e B -8£®J : (e))J : ^ 1 (e)d v h B =0,A=l,2,3. ' 

Remark 11. Compatibility conditions of the system of equations for V A delivers 
the constitutive conditions on the balance system 

Remark 12. Equation 

d v h°-J : - 1 (?)d v h(?) = 

can be considered as an abstract infinitesimal version of the Gibbs relation of the 
locally equilibrium thermodynamics, [llj . 

Remark 13. It is interesting to see if the set £ is connected and if not, do we have 
the same or different SBL from the same h A l 
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12. Cattaneo heat propagation system. 

In this section we determine the form of all supplementary balance laws for the 
Cattaneo heat propagation system (7.6). Instead of solving the LL-equations (4.5) 
directly, we will use the exterior systems generated by vertical differentials of the 
flux components F A . 

Consider the heat propagation model containing the temperature d and heat 
flux q as the independent dynamical fields. y° = $ , y A — q A , A = 1, 2, 3. 

Balance equations of this model have the form 



\d t (pe)+div(q) = 0, 
\d t (Tq) + VA($) = -q. 

Second equation can be rewritten in the conventional form 



(12.1) 



d t (Tq)+\-\7d = -q, 

where A = If coefficient A may depend on the density p, equation is more 
complex. 

Since p is not considered here as a dynamical variable, we merge it with the 
field e and from now on and till the end it will be omitted. On the other hand, in 
this model the the energy e depends on temperature d and on the heat flux q (see 
[13], Sec. 2. 1.2) or, by change of variables, temperature ■& = i?(e, q) will be considered 
as the function of dynamical variables. 

All variables may depend on x^ directly or though p. 

Cattaneo equation q + rdt(q) = — A ■ V?9 has the form of the vectorial balance 
law and, as a result there is no need for the constitutive relations to depend on the 
derivatives of the basic fields. No derivatives appears in the constitutive relation, 
therefore, this is the RET model. In the second equation there is a nonzero pro- 
duction H A = —q A . Model is homogeneous, there is no explicit dependence of any 
functions on t 1 x A . 

Constitutive relation specify dependence of the internal energy e on iJ,? and 
possible dependence of coefficients r, A on the temperature. Simplest case is the 
linear relation e = fci?, but for our purposes it is too restrictive, see [14], Sec. 2.1. 

Vertical variables are here $,q A - 

For the components of constitutive relation we have 

F°=e(#,q A ), F°=T($)q A , F B = q B , F§ = 8§A(#), A,B,C = 1, 2, 3. (12.2) 



We start with the i x p matrix 



Ff 



rq L 
A(0) 







rq 


A(0) 




rq" 



\q° A(i?)/ 

Assuming that coefficients r and the function A are independent on the vertical 
variables q A we get the vertical differentials 

v dd + rdq 1 Tf,q 2 dd + rdq 2 T{,q 2 dd + rdq 3 \ 
A#d$ 



d v F 



( e#dd + e q A 
dq 1 



V 



dq 2 
dq 3 



A 6 dtf 






A#d$ 



J 
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Left side of LL-subsystem for n = is 
H°(X) = (A , A 1 , A 2 , A 3 ) • 



/ £i5 






e 9 3^ 




fd a oK°\ 




T 










d a iK° 







T 







d a 2K° 














\d a3 K°J 



(12.3) 



As a result, LL subsystem for fj, = takes the form 

fA^ + TtfAV^ ,, 
{ A° V + A^r = 

For n = 1 we have in the same way 



, 4=1,2,3. 



(12.4) 



(A , A 1 , A 2 , A 3 )- 



Repeating this for A 
in the form 



(o 


1 





°\ 


A*, 























\o 








<V 



K 1 



'A^A 1 = if 1 ,,, 
= ^. 



(12.5) 







2, 3 we get the subsystems of LL-system for 4 = 1,2,3 



' Ai)X A 

X° = K / 
i 

= K A A , 



(12.6) 



) 4—1,2,3. 



Looking at systems (12.4-6) we see that if we make the change of variables: 
d = A(i?) then the systems equations we get takes the form ( wherever is the 
derivative by ■& we multiply this equation by A^) 



A e o A + X a t = K? A 




1,2,3. 



(12.7) 



Second subsystem is equivalent to the relation 

d v K A = X A dd + X°dq A . 
Integrability conditions of these form imply K A = K A ($,q A ) and 
K A A =X°, A = 1, 2, 3 => X° = A°(0). 
Integrating equation K A A = A (•&) by q® we get 

A^ A = A°(i9)g A + A >A (^). 
First equation of each system now takes the form 

X A =Kf = X° d q A + K A d (#). 

Substituting these expressions for X A into the 0-th system 



(12.8) 
(12.9) 



% A %a + X a t = K° qA 



1,2,3, 
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we get 

(K°, = X°er, + T-JXlWaW 2 + K A Jti)a A ). 

A = 1,2,3, (12.10) 



K s ^\°e d + T d (\l\\ q r + K^) q A ), 



where \\q\\ 2 = J2a1 A2 



a ' 

Integrating A-th equation by q A and comparing results for different A we obtain 



the following representation 



K° = A°e + r(d)[l\l\\q\\ 2 + K%d)q A ] + f{$) (12.11) 



for some function /($). 

Calculate derivative by d in the last formula for K° and subtract the first formula 
of the previous system. We get 

o = + r(m^°M 2 + k%^ a )]^ \\^hf + ^W- ( 12J2 ) 

This is the compatibility condition for the system (12.11) for K°. As such, it is 
realization of the general compatibility system (6.2). 

Take q A — in the last equation, i.e. consider the case where there are no heat 
flux. Then the internal energy reduces to its equilibrium value e^^d) and we get 
/#($) = —X°~e eq . Integrating here we find 



W) = fo- f X^{s)e e i{s)ds. (12.13) 
Substituting this value for / into the previous formula and we get expressions for 

A o e _ f X y qds + T( ^)[i A o || 5 ||2 + Rp)^ + /oj 

= X°(#)q A + K A (#), A = 1,2, 3. 

In addition to this, from (12.12) and obtained expression for /(i?), we get the 
expression for internal energy 



!_ „_„2 T@) 



e = <■**(#) + -T 9 \\q\\- 



2 A °(0) 



(12.15) 



This expression present the restriction to the constitutive relations in 
Cattaneo model placed on it by the entropy principle. 

Zero-th main field A is an arbitrary function of d while X A is given by (12.10): 

\ A = (\°y + Kffi). (12.16) 

Using this we find the source term 

Q = A A n A = —X A q A = -(A° |M| 2 + K A (d)q A ). (12.17) 

Now we combine obtained expressions for components of a secondary balance 
law. We have to take into account that (see (12.7)) that the LL-system of relations 
defines K» only mod C°°(X) (RET case!). This means first of all that all the 
functions may depend explicitly on x^ . For energy e, field A(i?) and the coefficient 
r this dependence is determined by constitutive relations and is, therefore, fixed. 
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Looking at (12.15) we see that the coefficients f linear and quadratic by q A are also 
defined by the constitutive relation, i.e. in the representation 



e = e«{&)+nW\\ q f+M A @)q A = e^)+ir J g || 2 --^ 
coefficients 



Kuh\\ 2 + KU^)q A 



1 

2 ,00 M * M 1 ",00 

(12.18) 



M,x) = -T S - -^L\%^ M A = (12.19) 
2 ^ 2A°-(i?) ^ A°-(i?) ' 

are defined by the CR - by expression of internal energy as the quadratic function 
of the heat flux. 

Rewriting the first relation we get 

{ln(Xl)) i9 = ln(r)j 2^5 => i„(A|) = in(r) +^-^f ^{s)ds => 

=> A° = are- 2 /* a = e fc ° > 0. (12.20) 

From this relation we find 

/0 
[re' 2 r ^) ds ]du (12.21) 

Using this in the second formula (12.19) we get the expression for coefficients K A 

K A = —M A ■ ^ = -M A ae-*f S => 
,00 T (tf) 

/$ r w 

dw / [M A {u)e- 2 S u ^^ ds ]du. 

(12.22) 

Functions K A (d) are defined by the second formula in the second line. 
Thus, functions \$,K A $ are defined by the constitutive relations while coeffi- 
cients a > 0, a , k A , m A are arbitrary functions of . 

Remark 14. It would be interesting to chose lower limit in the integrals in previous 
formulas from physical reasons. Then the coefficients a , k A , m A might have some 
physical meaning too. 

Combine obtained results and returning to the variable $ we get (here we re- 
peatedly use the relation = fi^f^ = O^) -1 /^ = A ^V,0) 

K° = \°e f \yUs + t$)[\\% |M| 2 + K A {#)q A ] + / = 
= (a + aA°)e -af \%e^ds + ^ [f A° ||g|| 2 + (k, d k A + aK A (#))q A ] + /„, 
K A = \°(d)q A + K A (d) = (a + aA°(?9)) g A + k A A($) + m A + aK A {{)), A = 1,2, 3. 
Q = -X A q A = -{\l\\q\\ 2 + kA{d)q A ) = -A^(Ag ||<z|| 2 + A^k A q A + aK A (#)q A ) = 
l-Aj(a\Uq\\ 2 + A,#k A q A + ak A (#)q A ). 

(12.23) 
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Collecting obtained results together we present obtained expressions for sec- 
ondary balance laws first in short form and then - in the form where relations 
(12.19-20) were used to separate the original balance laws from the general form 



K 1 

K 2 
K 3 



= a 



i 



w 



f\°e f \ ^ds + r(tf)A^[±A° || g || 2 + k A (d)q A ] + f \ 
X a ({))q 1 + K 1 ^) 
\ a {d)q 2 +K 2 (§) 
\ a {d)q 3 +K 3 {d) 
-A^(\°Jq\\ 2 + k A (#)q A ) 

( K%{°d)q A \ 

k 2 {§) 
+k 3 (§) 
\-k-;k A ^) q A ) \ 



-E fc " 



(r{d)q A \ 
8\A{d) 

6 A A(#) 



\ 



) 



X°(-d)q 1 
X°(tf)q 2 
\°(d)q 3 





ffo\ 




m 1 


+ 






m 3 




W 



Remark 15. Notice the correspondence between the tensor structure of the basic 
fields of Cattaneo system - one scalar field (temperature and one vector field 
(heat flux q A ,A = 1,2,3) and the structure of space SBC(C) of supplementary 
balance laws - elements of SBC{C) depend on one scalar function of temperature 
A°($) and one covector function of temperature k A - 



Returning to the variable ■& in the expression (12.15) and using the relation 



■d ~ a(i?) a dfl we get the expression for the internal energy 

„ 1,2 ^) 



2A 



o 



A, 



_A^=K-const £ eq^ _|_ < v ||g||2 _ 

2k 



II4II 



9/ 



Ah\\ 



k%q A 



(12.25) 



Notice that for A = 0, balance law given by the 4th column in (12.24) vanish. 
The same is true for deformations of the Cattaneo equation defined by the third 
column. 

First and second balance laws in the system (12.24) are the balance laws of the 
original Cattaneo system. Last one is the trivial balance law (see Sec. 3). Third 
column gives the balance law 



dt 



A°e- / A° 



\e e «ds + r(tf)A^[-AS||g|| 2 + K A (#)q A ] 



+8 x a 



\°(#)q A + k A (#) 



= -Aj(\°A4 2 + K A ^)Q- A )- (12-26) 
Source/production term in this equation has the form 
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- A J(X°J4 2 + K%{$)q A ) = -A-jA{h\\ 2 + 



AS 



1 A ) = 



- A^A 

— _A ,0 A # 



2AS 



(12.27) 



For a fixed i? this expression may have constant sign as the function of q A if and 
only if all K A ("&) — 0. Therefore this is possible only if the internal energy (12.25) 
has the form 

A, J 



2A 



2A° 



(12.28) 



with some function A°(i?), Cattaneo model has the supplementary balance law 



A°e 



A^e^ + lrWA^ 1 ^!!^! 2 



0* 



A (0)9' 



-A:/A°i| g | 



(12.29) 

with the production term that may have constant sign - nonnegative, 
provided 

A^Ag ^ 0. (12.30) 
This inequality is the II law of thermodynamics for Cattaneo heat prop- 
agation model. 

We collect obtained results in the following 

Theorem 7. (1) For the Cattaneo hear propagation balance system (12.1) com- 
patible with the entropy principle the internal energy has the form (12.25). 
All supplementary balance laws are listed in (12.24). New supplementary 
balance laws depend on the 4 functions of temperature - A° (i?) , K A (#) , A = 
1,2,3. 

(2) Additional balance law (12.26) given by the sum of third and forth columns 
in (12.24) has the nonnegative production term if and only if the internal 
energy e is given by the formula (12.28) and, in addition, the condition 
(12.30) holds. 

13. Conclusion. 

In this work we study the mathematical aspects of the development in the contin- 
uum thermodynamics that started with the pioneering works of B.Coleman, W.Noll 
and I. Muller in 60th of XX cent, and got its further development mostly in the 
works of G. Boillat, I-Shis Liu and T.Ruggeri. This development is known as the 
" Entropy Principle" . It combined in itself the structural requirement on the form 
of balance laws of the thermodynamical system (denote it (C)) and on the entropy 
balance law with the convexity condition of the entropy density. 

First of these requirements has pure mathematical form defining so called "sup- 
plementary balance laws" (SBL) associated with the original balance system. Space 
of SBL can be considered as a kind of natural " closure" of the original balance sys- 
tem. This vector space includes: original balance laws, the entropy balance, the 
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balance laws corresponding to the symmetries of the balance system and some other 
balance equations (see, for instance, Example 1). Thus, the space of supplementary 
balance laws carries important information about the thermodynamical system (C) 
and deserves the attention. 

In this text we study the case of Rational Extended Thermodynamics where 
densities, fluxes and production terms of the balance system do not depend on 
the derivatives of physical fields y % . We revisit the formalism of "main fields", 
Lagrange-Liu equations and dual formulation of the balance system in terms more 
formal then it is done in physical literature. Our main goal here was to obtain and 
start studying the defining system of equations for the density h° of a supplementary 
balance law. This overdetermined linear system of PDE of second order contains in 
itself equations determining all the densities h and with them, due to the formalism 
of RET, the fluxes and sources of SBL. In addition to this, this system contains, 
in the form of solvability conditions, the constitutive restrictions on the balance 
equations of the original balance system. 

We illustrate our results by some simple examples of balance system and by 
describing all the supplementary balance laws and the constitutive restrictions for 
the Cattaneo heat propagation system. 

Study of the space of supplementary balance laws for other physical systems (5 
fields fluid, Elasticity, 13 fields) is in the process and will be presented elsewhere. 

14. Appendix. 

14.1. Appendix I. Solvability of systems of PDE.. Here we remind the fol- 
lowing definition ([23], (2.70)). 

Definition 5. (1) A system of differential equations A(x /i , y l {x), z\\k ^ k) = 
0, is called locally solvable at a point zq = (xq, yo, . . . z\ ) G J k {ir) if 
there exists a smooth solution y — s{x) defined in a neighborhood of the 
point xo which has prescribed "initial condition" j k s(xo) = zq. 

(2) A system called non- degenerate at a point zq G Jp{^) if it is locally solvable 
at z Q and is of maximal rank at this point. 

(3) A system is called locally solvable (non- degenerate) in a domain D C X if 
it is locally solvable (non- degenerate) at each point of domain D. 

14.2. Appendix II. Symmetric hyperbolic systems. 
Definition 6. A system of linear equations 

^ = IT ^ 
for a vector function y{x^) is called symmetrical hyperbolic system if 

(1) Matrices A fl (x,y) are symmetrical, 

(2) Matrix A is positive definite. 

14.3. Appendix III. Cartan Lemma. 

Proposition 8. (Cartan Lemma. |33j . . Thm.4-4) Let V be a vector space over 
a field k. Let vectors Xi,i = l,...,p are linearly independent and for a vectors 
Vi G V,i = l,...,p, 

^x l Ay t = 0. 
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Then 

V 
3 = 1 
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